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Diagrammatic Monte Carlo approach is applied to a problem of a single spin-down fermion reso- 
nantly interacting (we consider the universal limit of short-range interaction potential with divergent 
s-wave scattering length) with the sea of ideal spin-up fermions. On one hand, we develop a generic, 
sign-problem tolerant, method of exact numerical solution of polaron-type models. On the other 
hand, our solution is important for understanding the phase diagram and properties of the BCS-BEC 
crossover in the strongly imbalanced regime. This is the first, and possibly characteristic, example of 
how the Monte Carlo approach can be applied to a divergent sign-alternating diagrammatic series. 
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In this Letter we introduce a novel technique for stud- 
ies of polaron type models which is equally important 
from both physical and technical points of view. On the 
physical side, we find a controUably-accurate numeric so- 
lution for the problem of a single spin-down fermion reso- 
nantly interacting with the sea of ideal spin-up fermions. 
This problem naturally arises in studies of the BCS-BEC 
crossover in the strongly imbalanced regime [l|, 01- In 
the case of stability of the dilute spin-down sub-system 
in the ideal spin-up fermi-sea, the solution of the single- 
particle problem would naturally define the phase dia- 
gram in the vicinity of the multicritical point (so-called 
M-point), discussed recently by Sachdev and Yang 3, 
where four different phases meat. Even if the M-point is 
thermodynamically unstable — as is actually suggested by 
the analysis of Giorgini and collaborators on the basis of 
their fixed-node Monte Carlo simulations [1] — the single- 
particle M-point, corresponding to the critical interaction 
strength at which the spin-down fermion forms a bound 
state with a spin-up fermion thus becoming a spin-zero 
composite boson, is of interest on its own. Also, for mak- 
ing a definitive theoretical conclusion about thermody- 
namic (in)stability of the dilute spin-down sub-system 
one still has to find an unbiased solution to the single- 
particle problem. 

On the technical side, we solve the Fermi-polaron 
problem by diagrammatic Monte Carlo (MC) technique, 
which proved very efficient for clcctron-phonon polaron 
problems Though the diagrammatic series for the 
Fermi-polaron are quite similar, the crucial difference is 
that in the Fermi system we have to deal with the sign- 
alternating, divergent (at least for strong coupling) series. 
Under these conditions a direct summation of all relevant 
Fcynman diagrams for the Green's functions is not possi- 
ble, and one has to develop additional tools for (i) reduc- 
ing the number of diagrams by calculating self-energies 
rather than Green's functions, and (ii) extrapolating MC 
results to the infinite diagram order for a divergent series. 
We find that the series for the Fermi-polaron are Riesz- 
summable and the numerically exact solution within the 
diagrammatic MC approach does exist. We believe that 
our findings are important in a much broader context 
since the diagrammatic MC approach to the many-body 



problem has essentially the same structure. 

In what follows, we will be using the term 'polaron' 
in a narrow sense, i.e. only for the unbound fermionic 
spin-down excitation. For the composite boson we will 
be using the term 'molecule'. If polaron is a well-defined 
elementary excitation, its energy E{p) can be extracted 
from the pole in the single-particle spin-down Green's 
function in the frequency-momentum representation, by 
solving the equation 

[Gf\uj = E{p), p)]-i = Si(c. = E{p), p) , (1) 

where G^^"^ is the vacuum Green's function for the spin- 
down particle and is the self-energy of spin-down par- 
ticle in the Fermi-sea of spin-up particles. In diagram- 
matic Monte Carlo we traditionally use the imaginary- 
time-momentum representation to circumvent the prob- 
lem of dealing with the singular structure of free propa- 
gators and alleviate the sign problem. Remarkably, in or- 
der to find the polaron energy we do not have to perform 
a numeric analytic continuation from imaginary to real 
frequencies, because in the imaginary-time-momentum 
representation Eq. ^ translates into 

/•OO 

[Gf\u; = E{p), p)]-i = / E^r, p)e^(P)-dr . (2) 

Jq 

The equivalence of Eqs. ([T]) and ^ readily follows from 
the fact that for the spin-down particle there is a freedom 
of choosing the 'chemical potential' (for a single par- 
ticle this is just a uniform external potential which does 
not affect its physical properties). This freedom can be 
utilized for fine-tuning fi^ = E{p), in which case the so- 
lution of Eq. ID) corresponds to zero frequency. Then, by 
noting that zero real frequency is the same as zero imag- 
inary frequency, and utilizing the trivial dependence of 
self-energy on fj^i, namely, Sx(t, p, /ij = ^^(t, p)e''i^, 
one proves Eq. 

We obtain self-energy (r, p) by means of diagram- 
matic MC which works with standard diagrammatic ex- 
pansions for Green's functions by interpreting them as 
a statistical ensemble. In our case, the diagrams are 
constructed from the following lines: (i) the already- 
mentioned vacuum propagator Gl"''(r, p) = e'^^i~^'^/^™)'^ 
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FIG. 1: The first-order diagram for E|, consisting of the 
T-matrix propagator (heavy dashed hne) and spin-up propa- 
gator (sohd arc). 




FIG. 2: The first diagram for molecular self-energy K, con- 
sisting of one T-matrix propagator (heavy dashed line), two 
spin-up propagators (solid arcs), and two spin-down propaga- 
tors (solid straight lines). 



where m is the particle mass; it will be represented by 
thin horizontal straight lines, (ii) the spin-up propaga- 
tors G|(r,p > kp) = e('=^-p'/2m)^ and G^{t,p < kp) = 
Q-{eF-p /2m)T (-^^jjQj-g a,nd ep are the Fermi momen- 
tum and Fermi energy of spin- up particles); they are de- 
picted with thin-line arcs, and (iii) the T-matrix prop- 
agator r(r, p). The specifics of the resonant problem is 
that T-matrix has to be considered as a separate dia- 
grammatic element which sums ladder diagrams for the 
short-range pair interaction potential. This summation 
takes the ultra-violet physics into account exactly and al- 
lows to express r(r, p) in terms of the s-scattering length 
a. The ladder structure of diagrams absorbed in the T- 
matrix explains why we treat it as a 'pair propagator' and 
depict it with a (dashed heavy) line. The exact expres- 
sion for r(r, p) in the zero-range limit in the frequency- 
momentum representation reads 



m 
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Ami] — n(p, 77) , (3) 
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As an illustration, in Fig. [T] we present the first-order 
diagram for E| . 

The energy of the molecule, i?m(p), cannot be ex- 
tracted from Eqs. ([1]), ([2]), because the corresponding 
pole exists only in the two-particle (spin- up + spin-down) 
or "pair" channel. Due to the fact that the resonant 
T-matrix is formally playing the role of a pair propaga- 
tor, the theory of the molecular pole is analogous to the 
Dyson-equation theory of the polaron pole. Diagrammat- 
ically, it arises from the summation of geometric series 
of products TKTKTK . . . where K is the T-irreducible 
diagram in the two-particle channel — a direct analog of 



self-energy of the single-particle channel (see Fig. [2] for 
an illustration). Correspondingly, the molecular energy 
is found by solving an equation 

/•oc 

T-\u; = i?™(p), p) = / K{t, p) e^-fP'-dr . (6) 

We found neither analytical expression for r(T, p), nor 
fast way of tabulating it with high accuracy using inverse 
Laplace transform of the frequency-momentum represen- 
tation Eq. ([3]). Instead, we applied recently-developed 
bold diagrammatic Monte Carlo (BMC) technique Q 
to obtain T[t,p) numerically, by relating it to the vac- 
uum T-matrix the analytic expression for which in (t, p)- 
representation is readily obtained. This was the first 
practical apphcation of BMC, see Ref. Q for details. 

The updates we use for sampling the diagrammatic se- 
ries do not differ dramatically from the ones described 
in literature, and we leave the corresponding discussion 
to a full-sized paper 0. Instead, we concentrate on 
the convergence issues. In contrast to previous exam- 
ples of diagrammatic Monte Carlo, the series for the 
resonant Fermi-polaron problem turns out to be diver- 
gent. However, the Cesaro-Riesz summation techniques 
fix the problem. For the quantity of interest — polaron 
or molecule self-energy — we construct the Cesaro-Riesz 



N=l ^N^N 



defined as sums 



partial sums, ^{N.^,) = ^ 
of all diagrams up to order TV* [the diagram order, A'^, is 
defined by the number of spin-up propagators] with the 
A^-th order terms being multiplied by the factor 
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(Cesaro-Riesz) . (7) 



Here 5 > is a fixed parameter {S ~ 1 corresponds to the 
Cesaro method). If the series is Riesz-summable, then 
the answer in the A* — > 00 limit does not depend on S. 
The freedom of choosing the value of Riesz's exponent S 
is used to optimize the convergence. 

For the A^*-truncated and reweighed series we first 
determine the polaron and molecule energies and then 
study their dependence on A^* as A^* 00. Figure[3]illus- 
trates the procedure. The odd/even oscillations are very 
pronounced for S ~ 1 hinting at the absence of conver- 
gence of the original series but are strongly suppressed for 
larger values of i5. With i5 = 4 we were not able to resolve 
odd-even oscillations, but the smoothness of the curve 
here comes at the expense of increased curvature, which 
renders the extrapolation to the 1/A^* limit vulner- 
able to systematic errors. Empirically, we constructed a 
sort of generalized Cesaro-Riesz factor F^*"* which leads 
to a much faster convergence (cf. Fig.[3]); in what follows 
we refer to it as an advanced summation technique 
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FIG. 3: The molecule energy (at kp a = \) as a function of 
truncation parameter A'^* for different summation techniques: 
Cesaro (open squares), Riesz 5 = 2 (filled circles, fitted with 
the parabola y = -2.6164 + 0.28013a; + 0.01638a;^), Riesz 
5 = 4 (open circles, fitted with the parabola y = —2.6190 + 
0.616352; — 0.3515a;^), advanced summation described in the 
text (stars). In the 5 = 2 case, the odd-even oscillations are 
strongly suppressed, as compared to Cesaro case, but still 
visible at higher resolution. 
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FIG. 4: Summing the Grandi series, 1-1 + 1-1 + . .. = 1/2, 
by Riesz method and by the advanced method described in 
the text. The solid lines are to guide an eye. 

For the molecular energy, the advanced summation 
demonstrates no visible linear slope in the l/A^* — > 
limit. For the polaron energy, there is a small linear 
slope, but without visible curvature, see Fig. [S] 

To ensure that both the code and data processing are 
error free, we check our answer for the molecule ground- 
state energy against the asymptotic solution 

= \ - £f + 7:S-"T (fc^a ^ 0) , (10) 

in which the first term in the r.h.s. is the binding energy 



FIG. 5: The polaron energy (at the unitarity point a^^ = 0) 
as a function of truncation parameter A'', for the advance 
summation technique. The asymptotic behavior at l/Af* — » 
is perfectly fit by a straight line. 



of molecule in vacuum, the second term reflects finite 
chemical potential for spin- up fermions, and the third 
term comes from the mean-field interaction between the 
compact molecule and the remaining spin- up particles. 
Correspondingly, a « 1.18a is the molecule-fermion s- 
scattering length 3. The value of a is based on the 
non-perturbative solution of the three-body problem, and 
thus provides a robust test to the entire numerical pro- 
cedure of sampling divergent sign-alternating diagram- 
matic series. Our data are in a perfect agreement with 
the a w 1.18a result within the statistical uncertainty of 
the order of 5%. 

In Fig. [6l we show polaron and molecule energies in 
the region of kpa ~ 1 where the nature of the quasi- 
particle state changes. The M-point is found to be 
at {kpa)c = 1.11(2). The values of both polaron and 
molecule energies are in excellent agreement with the 
fixed- note Monte Carlo simulations 01 • Interestingly, 
the polaron self-energy is nearly exhausted by the first- 
order diagram considered in Ref. 0, see also Fig. [51 

A comment is in order here on why the polaron 
(molecule) remains a good elementary excitation at 
kpa < {kpa)c (kpa > (kpa)c) and the M-point is an 
exact crossing. In the limit of kpa — > {kpa)c, this is 
guaranteed by the phase-space argument. The conserva- 
tion of energy, momentum, and spin projection dictate 
that the leading decay channel involves four quasiparti- 
cles in the final state: the polaron decays into molecule, 
two holes and one spin-up particle; the molecule decays 
into a polaron, one hole, and two spin-up particles. Cor- 
respondingly, the decay width — proportional to the four- 
particle phase-space volume allowed by the conservation 
laws — gets negligibly small as compared to the energy dif- 
ference \Em — Ep\ at kpa {kpa)c- Our numeric results 
show that \Ep — Em\ «C \Ep\, \Em\ even when kpa is sub- 
stantially different from (kpa)c- Since the phase-space 
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FIG. 6: Polaron (circles) and molecule (triangles) energies 
(in units of ef) as functions of kpa.. The dashed line on the 
lower panel corresponds to Eq. (fTO)l . 

volume is determined by \Ey — -Eml, this explains why 
the width of the decaying quasiparticle remains small and 
goes beyond our numeric resolution. 

The data for the effective mass is presented in Fig. (O . 
As expected, at the crossing point the effective mass 
curve is discontinuous. Note that good agreement with 
Eq. pn]) for Em at all couplings up to the Af-point is 
slightly misleading since Eq. pO|) assumes that molecules 



are compact and their mass is 2m. The actual effective 
mass is significantly enhanced in the vicinity of the M- 
point. 

Summarizing, the problem of resonant Fermi-polaron 
has been solved by diagrammatic MC. In particular, the 
point where the groundstate switches from the single- 
particle (fermionic) sector to two-particle (bosonic) sec- 
tor is found to be at kpa = 1.11(2). On the computa- 
tional side, it has been discovered that while the diagram- 
matic series is divergent, the fermionic sign-alternation of 
the diagrams renders the series summable by the Cesaro- 
Riesz-type methods and suitable for the diagrammatic 
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FIG. 7: Polaron (circles) and molecule (triangles) effective 
mass (in units of m) as functions of kpa. The vertical dot- 
ted line stands for (kFOjc = 1.11. The dashed line is the 
contribution from the first diagram [^. 

MC routine. Whether this sign blessing is a specific fea- 
ture of the problem solved here, or a generic feature of 
fermionic and, possibly, other diagrammatic expansions, 
is a major question which we will address in future. 
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